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;ii ' Abstract 



We motivate and analyze the weak-field limit of a non-analytical Lagrangian for the 
gravitational field. After investigating the parameter space of the model, we impose 



00 I constraints on the parameters characterizing this class of theories imposed by Solar- 



System data, i.e. we establish the validity range where this solution applies and refine 

t:j- ■ the constraints by the comparison with planetary orbits. As a result, we claim that this 

[^ . class of models is viable within different astrophysical scales. 
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1 Basic Statements 

Although General Relativity (GR) is a well settled-down theory for the description of 
geometrodynamics, nevertheless, over the last decades, a wide number of approaches 
have been developed to generalize it. These extended points of view are aimed not only 
at addressing the quantization of the gravitational field, but also at establishing proper 
deformations of the Einsteinian dynamics towards the space-time singularities. 
From the very beginning, the possibility to reformulate GR by using a generic function 
of the Ricci scalar (see, for example, [1] for a recent review and the references therein) 
has appeared as a natural issue offered by the fundamental principles established by 
Einstein. Indeed, this kind of extended Einstein-Hilbert (EH) Lagrangian preserves 
the fundamental features at the basis of General Relativity and Equivalence principles, 
and induces a direct modification of the dynamics only (at least, as far as the metric 
approach is concerned). Recently, modified fiTZ) gravity has acquired interest in view 
of the possibility to describe within this redefined dynamics some unexpected features 
observed on the large-scale structure of the Universe, i.e. the Pioneer anomaly, the 
galaxy rotation curve behaviour, the Universe acceleration, and, finally, the removal of 
singular behaviours of the gravitational field [2]. 

However, it is important to remark that any modification of the EH Lagrangian is 
refiected onto a deformed gravitational-field dynamics at any length scale investigated 
or observed. Thus, the success of such fiJZ) gravity in the solution of a specific problem 
has to match consistency with observation in other length scales [31 HI [5]. In this respect, 
we emphasize that the small value of the present curvature of the Universe [^ leads us 
to believe that, independently of its specific functional form, the fiJZ) term must be 
regarded as a lower-order expansion in the Ricci scalar. On the other hand, it is easily 
understood that the peculiarities of such an expansion will be extremely sensitive of the 
morphology of the deformed Lagrangian. 

The most immediate generalization is of course to deal with a function of the Ricci scalar 
analytical in the point i? = 0, so that its Taylor expansion holds [7j. This approach is 
equivalent to deal with a polynomial form [HI [9], whose free parameters are available to 
fit the observed phenomena on different sectors of investigation. Despite the appealing 
profile of such a choice, it is extremely important to observe that it is not the most 
general case, since real (non-integer) exponents of the Ricci scalar are in principle on 
the same footing as the simplest case [TOl [H]. In this letter, we will concentrate our 
attention on such an open issue, and we will develop a modified theory of the form 
fiJZ) =1Z + 'jTZ'^ , where 7 and P are two free parameters to be constrained at physical 
level [5l [3]. In particular, the details of our model will lead us to deal with rational non- 
integer numbers for /3, and to restrict it in the most appropriate interval for physical 
interpretation at low curvature, 2 < /3 < 3. 

We analyze the weak-field limit of our modified Lagrangian and derive the corresponding 
spherically-symmetric field equations. Here we retain only those non-integer powers in 
the Ricci scalar, which represent the dominant effect after the linear approximation. The 
choice 2 < /3 < 3 allows us to distinguish the dominant non-Einsteinian terms from the 
non-linear ones of GR. The explicit solution of the derived system is found, and its main 



features are analyzed in view of possible constraints from observational datafU [5], fT2] . 
The main issue of our treatment is to demonstrate that our fiTZ) model is appropriate 
both to fulfill Solar-System constraints and to provide a significant break-down of the 
Newton law at galactic scales. The paper is organized as follows. In section two, we 
provide the details of the model. In section three, we analyze the weak-field limit of the 
fiJZ) theory. Section four is eventually devoted to phenomenological estimations, which 
fix the constraints on the model. Discussion and conclusion follow. 

2 Non-analytical f{Tl) model 

We consider the following modified gravitational action in the Jordan frame 

s = - ife / d'xy^fin) = -^ J d'xy^ (7^ + 77^^) , (i) 

where /5 is dimensionless, and the parameter 7 has the dimensions of length 2/^-2 y\Jq 
can define the characteristic length scale of our model as L^ =| 7 j^/^^^^^) j^ jg straight- 
forward to verify that ([T]) is non-analytical in 7?. = for non-integer, rational /3, i.e. 
its Taylor expansion in the vicinity of 7^ = does not hold. Similar models have been 
investigated, for example, in [10], mostly from a cosmological point of view; in [H], the 
possibility of an irrational exponent is also envisaged. 

We can gain further information on the value of /3 by analyzing the conditions that al- 
low for a consistent weak-field stationary limit. Furthermore, the study of the potential 
of the scalar field in the Einstein frame will define the whole parameter space {P'j}- 
Having in mind to investigate the weak field limit of our theory to pursue its prediction 
at Solar-System scales, we can decompose the corresponding metric as g^jy = 7]^,^ + h^^, 
where /i^j, is a small (for our case, static) perturbation of the Minkowskian metric rj^^. 
In this limit, the Einstein equations read 

R^. - lv>.uR - rf3R^-';^, +^/3r]^,nR^-' = (2a) 

R = 3r/3ni?^-^ (2b) 

where ; denotes covariant differentiation, D = g^^VpVa, and, for the sake of compact- 
ness, in this treatment we retain the usual notation of the curvature scalar R and the 
Ricci tensor R^i, calculated up to the weak-field limit, as we are evaluating the 0{h) 
orders. 

The structure of such field equations lead us to focus our attention on the restricted 
region of the parameter space 2 < /3 < 3. This choice is enforced by the fulfillment of 
the conditions by which all other terms are negligible with respect to the linear and the 
lowest-order non-Einsteinian ones. 

The modified gravitational action ([T|) can be cast in an equivalent scalar-tensor model in 
the Einstein frame [UlIIS], i.e. a scalar field ip minimally coupled to gravity, by means of 
the conformal transformation g^i^ -^ ^"^Qiiv, where if = —In/'. In addition, the further 

transformation ip -^ kcj), where k = y^f^ accounts for the right dimensions of a scalar 



field. Within this framework, the potential V{ip) describing the dynamics of the scalar 
field reads 

n^) = e^^n-l + /3)(-l±^)^^7 (3) 

The appearance of a minimum in the potential is crucial, since the cosmological im- 
plementation of this picture into an isotropic and homogeneous Universe suggests one 
that such a minimum becomes, sooner or later, an attractive stable configuration for the 
system U. 

From the analysis of the scalar potential, we see that, for 2 < /3 < 3, there are two 
extremal points 

and V{ip) behaves likes an attractor only for the cases 

1. 7 > 0, /5 > 2, /5 = 2^^, m,n E Z: ipo is a minimum, ipp is a maximum; 

2. 7 < 0, /? > 2, /? = 2^^, m,n E Z: ip^ is a maximum, v?/3 is a minimum. 

3 Weak-field limit 

From the analysis of the form of ( 12all2bp in the Jordan frame, we learn that it is possible 

to find a post-Newtonian solution by solving the set of (1 2 all 2 bl) up to next-to-leading 

order in /i, i.e. up to (9(/i^~^), and neglecting the Oih?) contributions only for the cases 

2</3<3. 

These considerations motivate the choice we claimed above concerning the restriction of 

the parameter (3 G (2,3). The most general spherically-symmetric line element in the 

weak-field approximation reads 

ds^ = (1 + ^)df - (1 - ^)rfr2 - dfi2^ (4) 

where $ and \1/ are the two generalized gravitational potentials, and dVl? is the solid- 
angle element. Within this framework, the modified Einstein equations fl2all2bp rewrite 



Rtt -\R- -il3V^R^~^ = (5a) 

Err + \R- T(3R'^-\rr +^ ^V^ R^~^ = (5b) 

Ree + \r^R - r/?ri?^-\, +^pr^V^R^-^ = (5c) 

R = -3T(3V'^R'^-\ (5d) 



^ In fact the energy density of the scalar field is monotonically damped by the Universe expansion and 
the later stages of the Universe must be compatible with a weak perturbation of the GR scheme. As 
a matter of fact, a constant value of (/? corresponds to an Einsteinian dynamics and its extremal-point 
character ensures the nearby stability of this configuration. 



where 
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Here , denotes ordinary differentiation and V^ = ^ + ^^. 



System (I5atl5dl) is solved by 
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where the integration constant 6 has the dimensions of length, and the dimensionless 
integration constant a can be set equal to zero. The integration constant A has the 

2/3-2 

dimensions of length 2-/^ , and $^ and \E'/3 are dimensionless, accordingly. Moreover, we 
find that A is well-defined only in the case 7 < 0, /5 = 2n/{2m + l), i.e. only for configu- 
ration ([2]) of the scalar-tensor description. In this case, A> 0. It is remarkable that the 
limit 7 — > 0, apparently restoring GR [Uj, is mapped, in the solution, into the opposite 
case 7 — s> 00 |£|. Indeed, this situation in which the typical length scale associated to 
the parameter 7 is arbitrarily large would correspond to the physical intuition of an 
arbitrarily large scale for the relevance of the non-Einsteinian contributions. 



4 Solar-System constraints 

The most suitable arena where to evaluate the reliability and the validity range of the 
weak-field solution (I7atl7cp is, of course, the Solar System [HI [3]. For this reason, we 
can specify (JTbl) and (TTcll for the typical length scales involved in the problem. To this 
end, we split $ and "^ of (TTb]) and (TTcI) into two terms, the Newtonian part and the 
post-Newtonian one, respectively, i.e. $ = ^n + ^pn = rg/r + $^(r/L^)^^^~^^*^'^~^) and 



^ = \Ifjv + "^PN = —Tg/r + \Ef^(r/L^)^(^ ^^^^ ^), where the integration constant 6 in (JTbl 



Tel l results naturally as the Schwarzschild radius of the Sun (apart from the sign) r^ = 
2GMs/c^, Ms being the Solar mass. In fact, while the weak-field approximation of the 
Schwarzschild metric is valid within the range r^ « r < 00 because it is asymptotically 
flat, the post-Newtonian correction has the peculiar feature to diverge for r — > 00. It is 
therefore necessary to establish a validity range rmin « r « Vmax, where this solution 
is physically predictivepTS]. The definition of this validity range is strictly related to the 



parameters /? and L^ involved in the model. Stringent constraints for L^ can be found in 
the analysis of planetary orbital periods, which represent a severe test within the Solar 
System, because of the high precision of experimental data. 

4.1 Validity range 

Since we aim to provide a physical picture at least of the planetary region of the Solar 
System, we are led to require that ^pn and '^pn remain small perturbations with respect 
to $7v and '^ n, so that it is easy to recognize the absence of a minimal radius except for 
the condition r » r^. A maximum radius r^ax appears instead. In fact, the definition 
of r^ax requires further discussion. The typical distance r* corresponds to the request 

I <l>jv(r*) |~ $Piv(r*), I *jv(r*) |~ ^Pjv(r*). (8) 

For Tg « r « r*, the system obeys thus Newtonian physics, and experiences the 
post-Newtonian term as a correction. Another maximum distance r** can be defined, 
according to the request that the weak-field expansion in ^ should hold, regardless to 
the ratios ^pn/^n and '^pn/'^n- According to our definition of rmm, r** is defined as 

I $iv(r**) |« $pjv(r**) ~ 1, I ^iv(r**) |« ^pjv(r**) ~ 1. (9) 

We remark that that, within this scheme, r* and r** are defined as functions of /3 and 
Ly, i.e. r* = r*{(3,L^) and r** = r**{(3,L^) respectively 

/3-2 2/3-2 /3-2 

r* ~ {rj^p) 3^^ Lf -' , r** ~ L^/^f~^ ■ (10) 

It is important to stress that^ for the validity of our scheme, the condition r* » r^ 
must hold, i.e. L^ » ^f-^^^^^^~^\,. 

4.2 Planetary orbital periods 

At this level, neglecting the lower-order effects concerning the eccentricity of the plan- 
etary orbit, we can deal with the simple model of a planet moving on circular orbit 
around the Sun, and its orbital period T is given by T = 27r(r/a)^/^, a = ^d^/dr being 
the centripetal acceleration. For our model, from (ffbl) . we get 

^/^=(c]&7.-^/Ml + 2fz|$/3^^ . (11) 

Hence we evaluate the correction to the Keplerian period Tk = 27rr'^/^(G'Ms)^^/^, com- 
pare it with he experimental data of the period T^xp and its uncertainty 5Texp and then 
impose that the correction be smaller than the experimental uncertainty, i.e. 

3/3-4 

(12) 
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where r^ is the mean orbital distance of a given planet from the Sun. 



4.3 Numerical estimations 

We now specify the previous considerations for and intermediate value of /3, allowed by 
the analysis of the weak-field limit consistency and the scalar-tensor description. To 
provide proper numbers, we choose a typical (non-peculiar) value of the parameter /3, 
say /? = 8/3. 

High-precision measurements are nowadays available for the distances between Solar- 
System planets and the Sun, so that the relative error in the orbital period is extremely 
small. According to this fact, we specify our analysis for example for the Earth [3], for 
which Texp = 365.256363051dai/s, ST^^p = 5.0 ■ IQ-^^days and r^ = 149.6 ■ lO^km. For 
P = 8/3, (p]) implies L^ > 1.147466382 ■ lO^^Km. 

This lower bound for L^ can be plugged into (JHl) and ([9]) to get an estimation of the 
distances r* and r**. Since the constraints on the tt and the rr components of the 
metric tensor imply similar conditions, we restrict our analysis to the physically-relevant 
case of the tt component, and, after direct calculation, obtain r* ~ 1.6 ■ lO^^km and 
r** ~ 1.5 ■ W^km. 

Such values of r*, r** and L^ are essentially stable with respect to the range of | << 
j3 « 3. In fact, in the vicinity of the greatest value /3 = 3, we would obtain L^ > 
1.3-10^^/cm, r* ~ 4-10^°/cm, r** ~ 1.3-lO^^km. Our analysis clarifies how the predictions 
of the corresponding equations for the weak-field limit appear viable in view of the 
constraints arising from the Solar-System physics. Indeed, the lower bound for L^ does 
not represent a serious shortcoming of the model, as we are going to discuss. 

5 Discussion and Conclusions 

When adopting an f(jVj model to extend GR, we are addressing an appealing procedure 
to consistently derive unexpected features of standard geometrodynamics. However, the 
wide spectrum of possible effects allowed by such a reformulation could represent one of 
the shortcomings of the approach if no precise account of observational data is taken. In- 
deed, we have enough precision in the measurement of astrophysical systems to represent 
a severe test for the reliability of an extended approach for the gravitational-field dynam- 
ics. The present investigation has demonstrated how, in contrast with other analogous 
treatments [m [3], the non-analytical case 2 < /5 < 3 is not excluded by Solar-System 
tests. Even though we have restricted our analysis here to the uncertainty of the Earth 
period, nevertheless in sheds light on the viability of the theory unambiguously. The 
lower bound of L-y is almost compatible already with the Solar-System scale, since it 
would predict non-Newtonian effects for outer regions only. However, it is clear that 
the corresponding value of r**, where the theory would require a non-linear treatment, 
is manifestly incompatible with observations on the galactic scale. Thus, we can obtain 
a more reliable estimation for the fundamental scale L^ by requiring that the value of 
r* overlap the typical galactic scale on which the rotation curves manifest the fiatness 
of their behaviour, as outlined by 21cm Hydrogen line [Uj. Indeed, taking r* in corre- 
spondence of the length of about lOKpc ~ 3 ■ lO^^fcm and an internal mass (assumed 



as spherically distributed) of lO^^M^, we get L^ ~ 4 ■ lO^^fcm and r** ~ 5 ■ 10^^ km 
for P = 8/3. It is not among the purposes of this short communication to discuss the 
detailed predictions of our theories on the galactic scales (we will address the analysis 
of this question in further investigations) . What we wish to emphasize here is that our 
proposal for a suitable value of the characteristic length L^ can predict small changes 
on Solar-System observables, while it could drive relevant modifications in higher-scale 
astrophysics. 
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